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ABSTRACT 

We consider several types of boundary conditions in the context of time domain models for 
acoustic waves. Experiments with four different duct terminations (hardwall, free radiation, 
foam, wedge) were carried out in a wave duct from which reflection coefficients over a wide 
frequency range were measured. These reflection coefficients are used to estimate parameters 
in the time domain boundary conditions and a comparison of the relative merits of the models 
in describing the data is presented. Boundary conditions which yield a good fit of the model 
to the experimental data were found for all duct terminations except the wedge. 
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0091. This research was carried out while the first two authors were visiting scientists at the Institute for 
Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, 
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1. Introduction 

Traditional techniques for modeling the sound fields in ducts and enclosures utilize 
an impedance concept to model the boundary conditions. This approach is quite useful 
for harmonic sound fields as well as quantification of steady state, random sound fields 
in the frequency domain. However, in order to solve for a transient response, impedance 
concepts which are defined only under steady state conditions in the frequency domain 
have not proven useful. When useful, these impedance based models provide a way of 
quantifying the acoustic properties of the boundary as a function of the various param- 
eters characterizing the acoustic interface. Generally, an acoustic impedance boundary 
condition is found to be a function of frequency and the spatial position along the bound- 
ary. In addition, the boundary condition is based on assumptions of linearity as well as 
that the boundary surface is locally reacting [1, p. 257-270]. The first assumption is 
questionable at sufficiently high sound pressure levels and the latter applies for resonator 
designs but not for the more generally used bulk reacting materials such as foams and 
fiberglass. 

For time domain solutions, the imposition of such a boundary condition is awkward 
except under severe restrictions such as a hard wall or a completely absorptive boundary 
where the boundary is not a function of frequency. This is a common problem in deriving 
time domain numerical solutions (e.g., finite difference, finite element) for a general duct 
propagation problem [2]. Since a general time domain excitation may be characterized bj' 
a varied frequency content as well as spatial mode distribution, it is not straight forward 
to apply these impedance boundary conditions in time domain models. If a simple spatial 
structure such as an incident plane wave may be assumed, then it is known for a linear 
acoustic response that the Fourier transform of the reflection coefficient may be convolved 
with the incident wave to derive the reflected wave. This concept was used by Bolton [2] 
to quantify in a time domain measurement the impedance characteristics of foam type 
materials. 

In order to develop state- space models for acoustic control problems, it is necessary 
to include boundary conditions coupled directly with the equations of motion. This may 
be done, for example, by using simple continuum oscillator models to yield the frequency 
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domain characteristics of typical acoustic boundary conditions ([1, p.263-264], [3], [4]). It 
is the purpose of this work to examine the ability of such models to describe the acoustic 
boundary conditions for a free end, a hard termination, an absorptive termination and 
a typical foam acoustic treatment. This is done for a range of frequencies for which 
the incident wave field is plane but for which the impedance (and reflection coefficient) 
is not constant. Results are presented comparing the measured values as a function of 
frequency to those modeled using a best fit to the data. It is found that different models 
afford better representation of the data depending on the particular acoustic boundary 
condition. 

The efforts reported here are the first steps in the development of state space - time 
domain models for use in control design problems related to active control of noise in 
a closed cylinder (such as an aircraft fuselage). In such applications, one has (due to 
the small distances involved) negligible medium damping of the acoustic pressure fields. 
Since the major dissipative mechanism entails the partial absorption, partial reflection 
that occurs at the pressure field/wall interface, it is important in the control of the 
acoustic pressure field to model this dissipation accurately. 
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2. Boundary conditions and reflection coefficients. 

The physical quantities used for the description of acoustic wave motion in a fluid 
are the acoustic pressure p (being the deviation from the mean pressure at equilibrium) 
and the velocity potential d>. For waves of small amplitude, to first order both p and <f> 
satisfy the linear undamped wave equation with uniform speed of sound c in the fluid, 
where c 2 = IjpK , p being the equilibrium density of the fluid and K being the adiabatic 
compressibility of the fluid at equilibrium. The velocity potential is a complex valued 
function such that v(t,x ) = - V fhe fluid’s velocity at x at any time t. The 

acoustic pressure is related to the velocity potential by p(f,x) = p(f>t{t,x) [1, p. 243-257]. 

We consider three types of boundary conditions for the wave equation and derive 
the corresponding reflection coefficients for simple-harmonic waves in a one-dimensional 
wave-guide by insertion of the superposition of right and left propagating waves into the 
boundary conditions. The reflection coefficients are of the form R = (z — pc)/(z + pc) 
where z is the acoustic impedance of the boundary surface ([1], p. 259 ff.). 

2.1 Oscillating boundaries. In [3], Beale considers a bounded region Cl C R 3 where 
the interaction of the fluid with the boundary material at the (sufficiently smooth) bound- 
ary dCl is modeled by a continuum of damped harmonic oscillators. Within Cl, the velocity 
potential <f>{t,x) satisfies the wave equation 

(2.1) 4> u (t,x) = c 2 \/ 2 <f>{t,x), xefi. 

For every x € dCl , the normal displacement 6(t, x ) of the boundary into the fixed domain 
Cl is assumed to be independent of other parts of the boundary surface (a surface of 
local reaction . [1, p. 260]). With the effective mass m(x), the resistance d(x) and the 
stiffness k(x) per unit area of the boundary surface, the interaction of the boundary and 
the interior pressure is described by 

(2.2) m(x)<5 tt (/,x) + d{x)6 t {t,x) + k(x)6(t,x) = -p<f> t {t,x), x G dCl. 

In addition, from the assumption that the surface is impenetrable by the fluid it follows 
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that the velocity of the fluid and of the boundary coincide, i.e., 

( 2 - 3 ) 6,(t,x) = <f> n (t,x), xedQ 

where <p n = |£ denotes the outward (to the domain fi) normal derivative of <f>. 

The system (2.1) - (2.3) can be formulated as an abstract Cauchy problem for the 
state u = (<j>,<t> t ,6,6 t ) in the space H = /f J (ft) x Z, 2 (H) x L 7 (dQ) x L 7 (dQ). Here, as usual, 
L 2 (M) denotes the Lebesgue space of square integrable functions on M C R 3 and //’( fl) 
is the quotient of the Sobolev space H l (Q) over the set of constant functions ( H m (fl ) is 
the space of functions having distributional derivatives up to order m in L 7 ( fl)). With 
the norm | • | H representing the energy of system (2.1) - (2.3) at any fixed time f, which 
is given by 

M H = f(p\ V <i>\ 2 + pc~ 2 \<f>t\ 2 )dx + j (*|£| 2 + m|<$,| 2 )c?5, 

H is a Hilbert space. For u(t) € ^(^4), where 

D(A ) = {(u!,u 2 ,u 3 ,v 4 ) € HI V 2 ui € T 2 (H),u 2 6 H l (Q),u ln = u 4 on 50}, 

the system (2.1) - (2.3) can be written in the form j t u{t) = Au(t ) with the linear operator 
A defined by 

A(u\, u 2 , U3, U4) = (u 2 ,c ^ 2 Uj, u 4, — (pu 2 + ku 3 -f- du, j) /m). 

In [3] the well-posedness of system (2.1) - (2.3) is established by proving that A with 
domain D{A) is the generator of a Co contraction semigroup on H, provided that m(x) > 
0,Ar(ar) > 0 and d(x) > 0 for each x G 50. An analysis of the spectrum is given and, 
in particular, Beale s considerations reveal that one cannot obtain an exponential decay 
bound for the semigroup even in the case where d > 0. 

We consider the one-dimensional case 0 = (-^,0) where the boundary at 0 is 
characterized by three oscillator constants m,d,k £ R. The general solution to the 
wave equation in (—^,0) is the superposition of right and left propagating waves. If 
(f>(t,x) = F(t — x/c) + G(t + x/c) and S(t) are (sufficiently smooth) solutions to the 
system (2.1) - (2.3) on (— C, 0), they satisfy the condition 

(2-4) mfiu(t) + d6 t (t) -f k8(t) = —pd t (t, 0) 
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(2.5) 


6 t (t) = <f> r (t,0). 


Computing 6 and 6 tt by integration (resp. differentiation) of (2.5) we obtain 
<5(0 = J\f\t) - G'{T))dT = -(F(0 - G(t))jc + const. 

6«(t) = -(F"(t)-G"( t))/c 

where the constant is set to zero because for F(t) = G(t) = 0 the boundary should be 
at equilibrium <5(0 = 0. Inserting this into (2.4), we find 


(2.6) mG"(t) + {d + pc)G\t) + kG(i) = mF"{t) + [d - pc)F\t ) + kF(t) 


for the coupling of F and G by the oscillator at x = 0. 

Suppose that F, the wave incident on the boundary, is a simple harmonic of frequency 
u>/27r, i.e., F(t — x/c ) = exp{iu,’(f — x/c)}. Then the right side of (2.6) is a harmonic forcing 
function for the linear oscillator (2.6). It follows (neglecting possible transients due to 
initial conditions that we do not specify here, see e.g., [1, p. 45]) that the steady state of 
G is harmonic with the same frequency u>/27r, i.e., G(t + x/c) = R(u>) exp{fu;(f + x/c )} 
with some complex constant F(u;). Inserting F, G and their derivatives into (2.6) we then 
obtain 


(OSCIL) 


R M 


mw ! — iu(d — pc) — k 
mu t 1 — iw(d -f pc) — k 


This formula for the reflection coefficient coincides with (6.3.8) of [1], where it is derived 
from equating the acoustic impedance p/ y j nc id en t with the impedance z(co) = iu;m + 
d + k/iu of the oscillating boundary, which implies R(u >) = (z(u) — pc)j(z(u) + pc). 

To sum up, in the one-dimensional version of the model (2.1) - (2.3), a simple har- 
monic wave is reflected at the boundary by reversion of its direction of propagation, 
multiplication of its amplitude with |F(uj)| and a shift of its phase by argF(w). The 
magnitude of R{u>) is smaller than one if and only if d > 0. In this case the oscillating 
boundary absorbs acoustic energy. In case d = 0 energy is conserved. In state space 
(time domain) terminology, we have (see [3]) that for d > 0, A is dissipative and the 
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infinitesimal generator of a C 0 contraction semigroup (i.e., 5(<) = e At where |5(i)|w < 1) 
so that acoustic energy can be absorbed at the boundary while if d = 0, S(t) is a unitary 
group and no energy is absorbed at the boundary (energy is conserved). However, even 
in the case that d > 0, one does not have an exponential decay bound for solutions (i.e., 
there does not exist an A > 0 such that |S(f)|// < Me~ Xi ). 

2.2 Damped elastic boundaries. For d = 0, k = 0 the model (2.2), (2.3) together 
with p = p4> t results in the boundary condition 

m{x)p n {t,x) + pp{t,x) = 0, x e dQ, 

for the acoustic pressure p. This is called a Robin or elastic boundary condition. To 
include dissipation it is extended by adding a p t term, which gives systems that are 
studied, for example, in [5], [6], [7]. In [7] the wave equation for the acoustic pressure 

(2.7) p n (t,x) = c 2 V 2 xGfi 

with damped elastic boundary conditions 

(2.8) ap(t,x) + 0p t (t,x) + cp n (t,x) = 0, x e dQ. 

is the model that is used for investigating active noise control techniques. 

For Q > 0, (3 > 0 the linear operator A defined by 

A(f,g) = (g,c>v 2 f) 


with domain 


D(A) = {(/,<?) € H 2 (Q) x H\n)\af + 0g + cf n = 0 on dSl) 


is the generator of a Co contraction semigroup corresponding to problem (2.7) - (2.8) with 

the state (p,p t ) in the space H } (Cl)x L 2 (Cl) taken with the /f 1 (f2)xT 2 (fl)-equi valent norm 

| • \h corresponding to the inner product <(/i, £?i), (/ 2 , <72))* = 7 / /1/2 + /v/i V / 2 + 

an n 
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\ f <7i< 7 2 (see [7]). Note that the semigroup is not a contraction in the usual H 1 (ft) x L 2 (ft) 
n 

norm. 

To obtain reflection coefficients for harmonic waves in the one-dimensional case ft = 
(— £, 0) we insert the superposition p(t,x) = F(t — x/c) + G(t + x/c ) into the condition 
op(t, 0) + Ppt(t, 0) + cp*(f,0) = 0, and this yields 

(1 + + aG(t) = (1 - /3)F'(t) - aF(t). 


Thus, if the incident w r ave F(t — x/c) = exp{iw(f — x/c)} is a simple harmonic, so is 
the reflected wave G(t -f x/c) = R(lo) exp{ito(t + x/c)} (aside from transients) with the 
reflection coefficient 


(ELAST) 


= 7 


?w(l — /?) — a 


iu>(l -f /?) + a 

This equals (z{u>) — pc)/(z(u>) + pc) with z(u) = iu>pc/(a + iujfl). 

This shows that (ELAST) models a boundary surface with specific acoustic impedance 
[1, p. 261] given by ( = zj pc — iu/{ot -f iu>0). The magnitude of the reflection is smaller 
than 1 if and only if /? > 0. In the case = 0 the impedance is purely imaginary and 
the amplitude of the reflected wave is equal to the amplitude of the incident one. In 
the state-space or time domain formulation we have that for /? > 0 the reflecting surface 
produces a contraction semigroup if the norm | • | h for //^(ft) x L 2 (ft) as indicated above 
is chosen. In actual fact, if j3 > 0 the damping is sufficiently strong so that one has an 
exponential bound in the 7 i induced operator norm: |5(t)|n ^ M 0 e~ Mot for M 0 . Po > 0; 
see [7]. 


2.3 Frequency-independent boundaries. For m = 0 and k = 0, the model (2.2), 
(2.3) reduces to 

(2.9) — d(x)4> n (t,x) = p<f> t (t, x), x € 5ft. 

This equation states that the acoustic impedance p/vj nc ;aent at the boundary equals d(x). 
Equation (2.1) with boundary conditions of the form (2.9) is the subject of many pub- 
lications throughout the literature on the wave equation. For aspects concerning well- 
posedness, decay of solutions or control problems see for example [8], [9], [10]. The 
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report [11] summarizes some results on the one-dimensional case Q = (— f, 0) of (2.1) 
with boundary conditions (which we refer to as impedance boundary conditions) of the 
form 


( 2 . 10 ) 


<Ai(*,0) + (c^(i,0) = 0 


with complex specific acoustic impedance ( of the boundary surface. In case Re £ > 0 
the operator 

A(w, p) = c(p, w") 

with domain 


D{A) = {(«j,p)|(u7,p) 6 H 2 (~e, 0) x H'(-(,Q),p + (w n = 0 at - f,0} 

generates a Co contraction semigroup S(t) corresponding to the state (c<f>, d t ) in a space 
He- The Hilbert space He is the product of the quotient space H l (— (, 0) (i.e., the 
quotient of H l over the set of constant functions) with L 2 (— (, 0). The norm \(w,p)\ 2 He = 
\ w '\i 2 + \p\ 2 L i in He corresponds to the acoustic energy of the state. 

Substituting <f>(t,x) = F(t — x/c) + G(t 4- x/c ) into (2.10), we obtain C(t) = {(£ — 
1)/(C + l)}T'(t) + const. (The constant here is set to 0, because F = 0 should imply 
G = 0 (passive boundary)). Thus, the reflection coefficient for this model 

(IMPED) R{w) = R= 

does not depend on the shape of the reflected waves. Moreover, Re ( > 0 is equivalent 
to |/?| < 1, with |/?| < 1 if and only if Re £ > 0. Thus, no acoustic energy is absorbed 
at the boundary if £ is purely imaginary. If ( = 1, then R = 0 and we have a totally 
absorbing boundary. In the state space formulation (see [11]) we find that for Re ( > 0, 
A is dissipative and generates a contraction semigroup 5(<). Indeed, for ( ^ ±1, we 
can argue that A is in fact a spectral operator and S(t) can be expanded in a Riesz 
basis of eigenfunctions. For Re ( > 0, energy is absorbed at the boundary and we 
obtain an exponential bound for S(t). For £ = 1, the spectrum of A is empty and 
moreover, we can argue that we have a totally absorbing boundary. Thus, we find that 
the state space formulation and the usual frequency domain considerations lead to similar 
characterizations in terms of the parameter (. 


8 



3. Experiments and Data Analysis 

3.1 Experimental Procedures. The procedure used to determine the reflection coef- 
ficient of the various duct termination conditions is similar to that outlined in reference 
[12]. The model for planar wave propagation in a duct represents the pressure anywhere 
in the duct as defined by the following equation. 

(3.1) p(t,x) = A(u?){exp[i(u>* - **)] + R ( u ) exp[i(uJt + kx)]} 

In this expression, t and x represent the temporal and spatial variables, respectively, 
p is the measured acoustic pressure relative to a common phase reference, A(u ;) is the 
unknown incident wave amplitude and R{yj) is the unknown complex reflection coefficient. 
The wave number, w/c, is represented by k and u> is the angular frequency. 

By measuring the pressure, p(t, xj), at a number of axial locations, xj, an over deter- 
mined set of equations may be formed using Equation (3.1) and solved in a least squares 
sense. The variables solved for in this analysis are the positive and negative complex 
wave amplitudes represented by A{u>) and A(u>)/?(u;). A schematic of the test configura- 
tion used in three of the test cases is shown in Figure 1. Ten microphone locations were 
utilized at five different axial locations. By making measurements in pairs on directly 
opposite sides of the duct and adding the two complex measurements, the contribution 
due to the first higher order mode is summed out. The resulting measurements at the 
five axial locations noted in the figure are used in conjunction with Equation (3.1). 

The data was acquired in the frequency domain using pseudo-random excitation of 
the acoustic sources shown in the left part of Figure 1. This harmonic excitation allowed 
the wave field to be defined with two hertz resolution over a bandwidth to 650hz. By 
exciting the duct with symmetrical excitation on opposite sides of the duct, the dominant 
wave generated was planar even above the cut-on frequency of the first cross mode at 
325hz. This in combination with the microphone averaging technique outlined above 
allowed the frequency range to extend to the cut-on frequency of the 2nd higher order 
mode at 650hz. 

For this configuration, three termination conditions were investigated. The first is 
a near hardwall condition attained by terminating the duct with a 0.5 inch thick reinforced 
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V 


Side view 



Figure 1 

aluminum plate. For a true hardwall termination, the reflection coefficient would be 
expected to be invariant with frequency and be purely real with a value of 1.0. This 
panel has a fundamental resonance at about 250hz and another at 450hz and therefore 
appears somewhat soft at these frequencies where minor variations from the hardwall 
condition are observed in the data. The data in Figure 2 (in all figures, the data will be 
given by a solid line) reveals a reflection coefficient with real part varying from 1.0 at 
20hz to about 0.8 at 600hz, while the imaginary part is near zero with variations between 
-0.1 and 0.1 in the range 20hz to 600hz. 
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The second case was for free radiation from an open duct. This case displays a 
frequency dependence as a plane wave at low frequency would reflect from a pressure 
release boundary with a purely real reflection coefficient of -1.0. However, at the high 
frequency limit, all of the propagating energy wave would beam out of the end of the open 
duct and the reflection coefficient should decrease to zero. The real part of the reflection 
coefficient data shown in Figure 3 exhibits this general trend. However, the imaginary 
part of the measured data varies over a positive range from 0.2 to 0.6. This behavior 
may be attributed to the presence of several reflecting surfaces outside of the duct but 
in the general vicinity of the duct exit. Also, the room was generally reverberant and 
may have exhibited some modal response. 

The final case tested in the duct of Figure 1 was an eggerate foam referred to as a 
wedge. This foam was 5.0Scm thick from tip to base and was backed by an additional 
10.16cm of closed cell foam. The back side of this foam was left open to the laboratory 
space. It is hard to anticipate the exact behavior of this type of termination condition. 
However, it is expected to have a generally complex reflection with a reflection coefficient 
that is close to 1 at low frequencies and that decreases as the frequency is increased. The 
data for this case is shown in Figure 4 and may be generally regarded to exhibit the 
correct trends. There is however a noticeable dip in the real part at 77hz. The cause for 
this low frequency behavior is not clear. 

In addition to the above cases, the material properties of a one inch thick acoustic 
foam backed by a hard surface were investigated using the impedance tube facility of 
reference [12]. This facility is designed for this type of measurement and better estimates 
of the material acoustic properties are to be expected. Data (depicted in Figure 5) from 
these experiments is referred to as the foam termination data below, and the general 
trends attributed to the wedge case above may also be attributed to this case. 

3.2 Computational Procedures. In order to evaluate the correspondence of the 
mathematical models of Section 2 with the measured data described in Section 3.1, 
the “difference” between the measurements for varying frequency / = w/2x and the 
models’ reflection coefficients R(co) was considered. By “difference” any kind of numerical 
variation between the measurements and the values R{uij) could be taken. For example, 
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if the phase of the reflected waves is of main interest (which in some sense is the case for 
noise suppression based on destructive interference), then “difference” could be based on 
deviation of the model’s phase shifts arg R from the measured ones. 

Here, we choose as “difference” the distance of the measurements to the models’ 
reflection coefficients as numbers in the complex plane. This seems to be a neutral choice 
in terms of a general comparison. To determine the best possible fit of the functions R(u>) 
to the data Rj, we minimized the functional 

SSQ = £ | Rj- *( Wj -)| 3 

3 - 1 

by variation of the parameters in R(l j). Here n is the number of measurements Rj at 
frequencies fj = u>j/ 2tt within the range of frequency considered. Because the data for 
the foam termination was taken in the facility of [12] at only discrete frequencies, there 
were only n = 23 frequency data points available between 20hz and lOOOhz. For the 
other three experiments the range of frequency considered is 20hz < fj < 600hz with 
n = 291. For each combination of a model with a duct termination, the minimal value of 
SSQ is an evaluation of the quality of the correspondence of the model with the data. 

Note that the data, that is derived from pressure measurements, is also applicable to 
the reflection coefficients for the velocity potential, because the right and left propagating 
components of the pressure waves are proportional to the time derivative of the corre- 
sponding components of the velocity potential i.e., if 4>{t,x) = F(f-x/c)+G(f+x/c) is the 
velocity potential, then p(t, x) = x) - p[F'{t - x/c) + G'{t + x/c)], which is consis- 
tent with Equation (3.1) if we choose F(£) = exp(iu>£) and (7(£) = exp(u<;£). 
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4. Results and Discussion 

The purpose of this work was to evaluate the ability of the proposed models to 
characterize a variety of real acoustic boundary conditions. As such there was no attempt 
to “smooth” the data by using analytical models or “perfect conditions” such as absolute 
hardwall or resonator impedance conditions in place of measured data. The data was 
collected for some general conditions which might be expected to be encountered in 
practice. 

For the minimization of the sum of squares SSQ over the parameters ( m,d,k ) E 
R 3 and (aj) E R 2 we used LMDIFl of the FORTRAN package MINPACK in the 
public-domain library NETLIB at Argonne National Laboratory. This routine is an 
implementation of the Levenberg-Marquardt algorithm with implicit scaling and optimal 
choice for the correction steps [13]. Starting from an initial estimate for the parameters, 
the algorithm converges to a local minimum. A systematic search with a variety of 
initial estimates and the use of graphics indicates that in all examples the numbers 
listed in Table 2 with corresponding parameters in Table 1 are the unique global minima 
of SSQ. In one case, for (OSCIL) and the hardwall data, a second local minimum 
10 3 x SSQ/n = 4.5075 is located at m = — 1708.4, d = 9.3346, k = —15384. 

Considering (IMPED), a simple algebraic argument shows that the gradient of SSQ 
with respect to the two real parameters (Re£, 7m£) E R 2 vanishes if and only if 
Re R = ]T)R e Rj/n and ImR = £ IrnRj/n. Thus, to get the optimal least-squares 
fit for (IMPED), we compute the mean values of Re Rj and ImRj and then £ according 

to C = (1 + R)/( 1 -R). 




Hard 

Free 

Wedge 

Foam 

Units 


m 

1028.0 

114.92 

-36.493 

85.621 

10 -3 kg/m 2 

OSCIL 

d 

8.6669 

0.1311 

0.2189 

1.4551 

10 3 kg/m 2 s 


k 

13607. 

-8.4717 

-788.62 

3851.9 

10 3 kg/m 2 s 2 

ELAST 

a 

-2.7985 

2178.4 

154.17 

-16.054 

1/s 


fi 

34.610 

761.85 

886.44 

146.15 

10" 3 

IMPED 

ReC, 

28.340 

0.3781 

0.8073 

4.8860 

1 


Im( 

-3978.6 

449.95 

717.72 

-3349.2 

10" 3 


Table 1: Minimizing parameters 
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The corresponding minimal sums of squares divided by the number of data n are given 
in Table 2. 


10 3 x SSQ/n 

Hard 

Free 

Wedge 

Foam 

OSCIL 

3.8983 

103.735 

51.586 

5.116 

ELAST 

4.8153 

7.819 

261.820 

54.729 

IMPED 

4.7972 

134.387 

137.038 

36.593 


Table 2: Residual sum of squares 

The experimental data (measured reflection coefficients * real and imaginary parts - 
as a function of frequency) are represented by solid lines in Figures 2 - 5 for hardwall, 
free radiation, wedge and foam terminations, respectively. The reflection coefficients 
corresponding to the models with boundary conditions (OSCIL), (ELAST), (IMPED) 
evaluated at the optimal parameters given in Table 1 are plotted in each figure by dashed, 
dotted, dashed-dotted lines, respectively. 
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Frequency in hz 
Figure 5b 
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As might be expected, none of the above three boundary conditions can be adjusted 
to all four of the boundary surfaces considered. The models for the impedance or re- 
flecting conditions are also limited in the degrees of freedom that they may exhibit. 
The impedance condition must by it’s definition be frequency independent. The elastic 
condition models only the prescribed frequency variation on the imaginary component. 
The oscillating boundary may model a somewhat more general condition as the resonant 
behavior may be tuned to any frequency and damped as necessary. However, one must 
keep in mind that there are only 3 real constants to adjust to model a general frequency 
dependence. It is clearly too much to expect that one single boundary condition will 
model any general acoustic boundary termination. 

But there are cases where the model curves fit to the data quite well: (OSCIL) for 
hardwall (/ < 450hz) and foam termination and (BLAST) for free radiation. 

The frequency-dependence of the measured reflection coefficients being a structural 
feature of the data (except for the hardwall with / < 450hz), the frequency-independent 
boundary condition (2.10) is not an appropriate model for the terminations free, wedge, 
and foam, unless only a narrow range of frequency is considered. 

In four examples the unconstrained minimization of SSQ renders negative parame- 
ters (m,fc,a < 0) which contradicts their physical definition and the assumptions used 
in the proofs of well-posedness of the models. Moreover, in these four cases the fit of the 
model functions to the data is not satisfying and would be worse if the positivity condi- 
tions were enforced in a minimization procedure with constraints. Thus, the oscillator 
boundary condition appears not applicable for free radiation or wedge termination nor 
is the damped elastic boundary condition suited for hardwall or foam termination. 

The set of data from the wedge termination experiment displays a pronounced mini- 
mum at 65hz and then decreases nearly linearly from about 150hz thru 600hz. The phase 
of the data set suffers the same anomaly near 65hz and then increases linearly from about 
40 degrees to 180 degrees at 600hz. Considering this variation, it is not surprising that 
the algebraically simple functions R(u) cannot mimic this data curve when only two or 
three parameters can be adjusted. But this is not the main reason for the difficulties 
with the wedge data. When consideration of the data is restricted to the range 250hz 
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< / < 600hz, where it is simply decreasing, the optimal mass and stiffness in (OSCIL) 
are again negative (the fit of Re R and Im R increases, 10 3 x SSQ/n = 7.049) and 
the optimal fit of (ELAST) again is not physically realistic (alpha, beta again positive, 
10 3 x SSQ/n = 76.009). This data was modeled assuming that the reflection plane was 
at the duct outlet which corresponds to the back of the foam layer. From an acoustic 
standpoint, the face of the wedges might be considered a more physically intuitive posi- 
tion to take as a reference. This however, only changes the phase response, increasing the 
slope significantly. Modeling the data in this way resulted in similar non-physical models 
and the data fit was even worse. These results suggest that the algebraic structure of 
the models considered here does not well represent the physical mechanisms inherent in 
the wedge boundary termination. 
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5. Conclusions 


In summary, the unconstrained least-squares fits suggest that among the boundary 
conditions considered, several might be appropriate for use in describing the reflection 
of harmonic waves by the duct terminations over the range of frequency considered. For 
the hardwall termination, the (OSCIL) boundary condition with physically reasonable 
parameter values provides a good agreement of reflection coefficients in the range / < 
450hz. The (OSCIL) boundary condition is also a reasonable choice for use with the foam 
termination. For the free radiation termination, the damped elastic conditions (ELAST) 
offer a good approximation of model to the experimental data. When the reflection of 
waves of arbitrary shape is to be modeled, a more detailed investigation of experimental 
data for such waves could give additional information on the quality of these boundary 
conditions over a wider range of frequency. 

For the wedge termination case the poor least-squares fits discussed above suggest 
that a more specific model for the interaction of the (non-flat) boundary surface with 
the interior field is needed to cover the experimentally observed phenomena. 
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